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ABSTRACT 

We show that gamma ray hne emission at ~0.4 MeV and ^0.2 MeV can be produced 
by Compton scattering of beamed radiation in the jets of Galactic black hole candidates. 
This mechanism has the novel feature of not invoking the presence of e~^-e~ pairs. To 
produce the two lines we employ a symmetric double sided jet with bulk flow velocity of 
about 0.5 c and incident beam radiation with a hard energy spectrum. We show that the 
two lines can be seen at viewing angle cosines relative to the jet ranging from 0.2 to 0.6. 
This comprises 40% of the total solid angle. In addition, the line radiation is approximately 
10% polarized. Depending on the bulk flow and viewing angle the model can produce lines 
at other energies as well. In particular a broad feature near 1 MeV can be seen by viewing 
the jet close to its axis. Our model can also accommodate single line spectra if the beamed 
gamma ray emission or the jets themselves are asymmetric. 

Subject headings: gamma rays: theory — radiation mechanisms: nonthermal — 
black hole: physics — line: formation 
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1. INTRODUCTION 

Transient gamma ray line emission at ~0.4 MeV was observed with the imaging gam- 
ma ray spectrometer SIGMA from the source 1E1740. 7-2942 located at about 1° from the 
Galactic center (Bouchet et al. 1991; Sunyaev et al. 1991) and Nova Muscae (Goldwurm et 
al. 1992; Sunyaev et al. 1992). This line was also seen with HEAO-l from an unidentified 
source 12° away from the Galactic center (Briggs 1991). The SIGMA observations of Nova 
Muscae revealed an additional line, at ^0.2 MeV which was also seen from the direction 
of the Galactic center (Leventhal &: MacCallum 1980). More recently, an observation the 
Galactic center region with HEXAGONE (a balloon borne Ge detector) revealed both 
the ~0.4 MeV and <0.2 MeV lines (Smith et al. 1993). The 1E1740. 7-2942 source could 
thus occasionally produce both lines, but because HEXAGONE was non-imaging, it is 
possible that the two lines were produced by different sources. In the same balloon flight 
HEXAGONE also detected the narrow 0.511 MeV line (Smith et al. 1993). This line was 
observed on many occasions from the direction of the Galactic center with other Ge detec- 
tors (Leventhal et al. 1978; 1980; 1993; Gehrels et al. 1991; Riegler et al. 1981; Albernhe 
et al. 1981) and with the OSSE instrument on CGRO (Purcell et al. 1993). Because of 
its precisely determined centroid, this line is undoubtedly due to e~^-e~ annihilation. 

The ~0.4 MeV line has also been interpreted as e~^-e~ annihilation radiation, but in 
this case the line center energy requires that the line be redshifted. The required shift could 
be gravitational, in which case the line must be formed within just a few Schwarzschild 
radii of a compact object, presumably a black hole. This implies that pair production 
and annihilation occur at essentially the same physical site, which leads to a problem 
because the temperature required to produce the pairs greatly exceeds the upper limit 
on the temperature set by the width of the ~0.4 MeV line. Moreover, for both 1E1740.7- 
2942 and Nova Muscae, the MeV photons which are required to produce the pairs have not 
been observed. The ^0.2 MeV line has been interpreted as Compton backscattered ~0.5 
MeV line emission from the inner edge of an accretion disk (Lingenfelter &: Hua 1991; Hua 
&: Lingenfelter 1993). This model produces a backscattered line of much smaller intensity 
than that of the primary line. Because, when observed simultaneously, the lines have 
roughly equal intensities, the model requires that the primary line be strongly attenuated. 
The attenuation could be due to the accretion disk itself. However, the upper limit on the 
temperature of the scattering medium [kT ^ 5 keV) set by the width of the ^0.2 MeV 
line is significantly lower than the temperature of the inner edge of the accretion disk. 
The redshift of the ~0.4 MeV line could also be due to the motion of the annihilating 
region, as in the model of Misra &: Melia (1993) where the pairs annihilate at the base 
of a jet emanating from the black hole. This model has the advantage of separating 
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the annihilation region from the pair production region, however it does not provide an 
explanation for the ^0.2 MeV line. 

In this paper we show that both the ~0.4 MeV and ^0.2 MeV lines could result from 
Compton scattering of coUimated hard photons by streaming plasma in a double sided 
jet. Radio observations have in fact revealed the presence of such a jet emanating from 
the 1E1740. 7-2942 source (Mirabel et al. 1991, 1992). The existence of coUimated high 
energy radiation is now well established from radio loud AGN (Fichtel et al. 1993) and the 
scattering of jet radiation may be responsible for sharp high energy turnovers in misaligned 
jet sources such as Centaurus A (Skibo, Dermer &: Kinzer 1994). In the scenario that we 
develop here transient gamma ray line emission is the result of an episode of enhanced 
bipolar energy release from a black hole in the form of plasma and hard photons coUimated 
into two oppositely directed jets. The optically thick material at the base of the jets scatters 
some of the beamed radiation into our line of sight. If the spectrum of the beamed radiation 
is sufficiently hard then the scattered radiation will exhibit spectral lines. The energies 
of the centroids of the two features determine both the orientation of the system with 
respect to the observer and the velocity of the scattering medium. The model can produce 
two lines with essentially equal strengths, it provides a relatively cool scattering medium 
consistent with the widths of the lines, it does not require pair production, and it does 
not predict the strong MeV emission which should be present if the pairs are produced 
isotropically by 7-7 collisions. The model predicts that both lines should be polarized. 
Even though the expected degree of polarization is quite modest, polarization data could 
be used to discriminate between our model and the models mentioned above which predict 
that the lines should be essentially unpolarized. Our model also predicts significant pair 
production. But these pairs are expected to annihilate outside the jets thereby contributing 
to only the narrow 0.511 MeV line emission. 



We consider the Compton scattering of a perfectly coUimated beam of photons with a 
power law spectrum directed in the positive z direction off a cold [kT ^5 keV) scattering 
medium moving with velocity /3c also along the positive z direction. In the optically thin 
regime the scattering rate per unit solid angle per unit energy in the observer's frame is 
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(see Skibo et al. 1994) where e is the energy of the scattered photons, /x = cos 6 is the cosine 
of the scattering angle, 7 = (1 — /3^)~2 ^ is the column density of the cold electrons 
along the z axis, Vf, = ^ lm^(? = 2.82 X 10~^^ cm and a and k are the spectral index 
and normalization factor, respectively, of the incident photon power law. The maximum 
energy of a scattered photon, obtained when the initial photon energy tends to infinity, is 

= 7(l+/5)(l-/x)- 

If a < 1 the spectrum given by equation (1) diverges as e approaches from below, 
resulting in a line-like feature at e*. This feature is caused by the pile-up just below of 
the large number of high energy photons implied by the hard incident spectrum. 

100 L ' J 



a 
5^ 




1 I , , . , I , , I I 

0.1 1.0 
e (MeV) 

Fig. 1. — The spectrum obtained by evaluating the sum of equation (1) and its 
reflection (/X —fl) for 7 = 1.3, a = and fl = 0.40 

If we consider two identical beams, one directed away from the origin in the positive 
z direction and the other in the negative z direction, then two lines result. The emergent 
spectrum is now the sum of equation (1) and its reflection /x — /x. In this case, the 
energies of the two lines are uniquely determined from the bulk velocity and the viewing 
angle. Conversely, given the energies of any two lines, the velocity and viewing angle are 
also uniquely determined. For example, if the two lines are at 0.40 and 0.17 MeV, then 
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from equation (2) 7 = 1.3 and /x = 0.40. In Fig. 1 we show the spectrum given by equation 
(1) and its reflection for this choice of 7 and /x, and with a = 0. The resultant hues are 
asymmetric, exhibiting a low energy excess and a precipitous cut off at e*(/x) and e*(— /x). 
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Fig. 2. — Spectra for various viewing angle cosines obtained from a Monte-Carlo 
simulation of photon scattering in a double-sided cylindrical jet of optical depth 
T|| =5 along the symmetry axis and T_|_ = 0.5 along its radius and moving with 
bulk Lorentz factor 7 = 1.2. To obtain the correct ordinates the ten spectra, from 
top to bottom, have to be multiplied by lO" to 10^, respectively. 

To calculate the spectrum and polarization of the scattered radiation for the case of a 
medium of flnite optical depth we resort to a Monte-Carlo approach. We take the scattering 
medium to be a cylinder of Thomson depths ry = (TtN^ along the z-axis and t±_ along its 
radius. In the simulation each photon is injected in the positive z direction at the base of the 
scattering cylinder with initial energy drawn at random from the energy distribution of the 
incident power law. To determine the scattering angle we numerically invert the cumulative 
probability distributions obtained from integrating the Klein- Nishina cross section over 
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scattering angle for fixed incident photon energies. The positions, directions and energies 
of each photon are tracked until they escape from the cylinder. We have for now ignored 
7-7 interactions which will remove some of the photons and generate pairs. The emergent 
photons are then binned in energy and angle cosine. To obtain photon fluxes with the 
motion of the scattering medium properly taken into account we apply the additional 
factor (1 — /3)/(l — /3/x). Provided the two jets are identical, each angle cosine bin contains 
photons scattered through an angle 6 from one jet and 180° — 6 from the other jet. 

The combined results for the two jets are shown in Fig. 2 for a = 0, maximum photon 
energy of 10 MeV, ry = 5, t± = 0.5 and 7 = 1.2. The photon fluxes at Earth from a source 
at a distance R are obtained by multiplying these spectra by R~'^ dNo / dt, where dNo/dt is 
the rate, measured in the observer's frame, at which photons are injected into each jet. As 
indicated above, 7 ~ 1.2 combined with a hard photon spectrum (e.g. a = 0) will produce 
line features at ~0.4 and ^0.2 MeV. Indeed, for 0.4 < /x < 0.5, the peaks of the two 
lines are at 0.16 and 0.43 MeV. The line center energies do depend on the observing angle. 
However, for /x in the range 0.2 to 0.6 (which comprises 40% of the total solid angle), the 
line energies are confined to the ranges 0.3-0.5 MeV and 0.15-0.18 MeV. Close to the z 
axis the spectrum is nearly identical to the original assumed flat spectrum and for small 
angles [6 ^45°) the spectrum exhibits a bump at ~1 MeV similar to that observed from 
Cygnus X-1 (Ling et al. 1987) and the Galactic center (Riegler et al. 1985). As 6 tends to 
90° the two lines merge into a single feature. 

We have calculated the polarization taking advantage of the fact that ry ^ 1 and 
T± ^1. For the assumed opacities, 64% of the photons scatter and 74% of these escape 
without further scattering. We then assign the polarization sin^ 6/{e/eo + eo/e — sin^ 6) to 
each singly scattered photon and zero polarization to multiply scattered photons. Here cq 
and e are the initial and final photon energies, and 6 is the scattering angle all evaluated 
in the rest frame of the scattering medium. In doing this we interpret each 'photon' of 
the simulation, not as a single quantum mechanical entity, but as an ensemble average. 
The calculated polarization is shown in Fig. 3 as a function of photon energy for the same 
ranges of observing angle cosines as given in Fig. 2. Even though polarizations as high 
as 35% are evident, maximal polarizations occur at energies outside the lines. The lines 
are formed predominantly by the downscatter of high energy photons, hence at the line 
energies the polarization is smaller. For example, for the 0.4 < /x < 0.5 spectrum, the 
polarization of both the lines is approximately 10%. 
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3. DISCUSSION 

To calculate the gamma ray luminosity L~f in the beam we rely on the ~0.4 MeV line 
observations which indicate fluxes of (0.2-1) xlO~^ photons s~^ cm~^ (Bouchet et al. 1991; 
Smith et al. 1993). Using the results of Fig. 2 with 0.4 < /x < 0.5, assuming a distance of 
8 kpc to the source and taking the line flux equal to 10~^ photons s 



^ cm ^, we obtain 



T ~ 



10 erg s . This equals the Eddington luminosity of a 10 object. 




£ (MeV) 



Fig. 3. — The degree of polarization of the scattered radiation. The individual panels 
correspond to the spectra in Figure 2. 

To maintain a stable scattering environment we require that the power in the plasma 
outflow Lp = n(3cA['y — l)mc^ exceed that in beamed radiation. Here n is the density, A 
is the cross sectional area of the outflow, and m is the proton mass in the case of ordinary 
matter or electron mass in the case of a pair plasma. Below we assume that Lp lO"^" erg 
s~^, which exceeds the Eddington luminosity for a 10 object. This should not present 
a problem because the radiation is beamed and the mass of the object may exceed 10 M©. 
We flrst consider the case where the scattering medium consists of ordinary matter. The 
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requirement that the Thomson opacities of the scattering cyhnder be large along its z axis, 
T|| > 1, and small along its radius rj^ < 1 leads to 

where L^dd is the Eddington luminosity, Rs is the Schwarzschild radius and ^ = t^\/t±. 
The right side inequality ensures that the cylinder is optically thick along the z axis while 
the left side keeps the cylinder optically thin along its radius. Using 7 = 1.2 and the above 
value of Lp, the right side limits the radius of the jet to ^140^Rs. In order not to have 
an excessively elongated scattering region, we take ^ < 10 which implies a radius less than 
1400 Rg. This is about 4 x 10^ cm for a 10 object. The left side inequality implies 
that the radius is larger than lAORg or 4 X 10^ cm because otherwise the transverse opacity 
will be too large. An intermediate value for the radius is then about 2 X 10^ cm, with a 
corresponding length about a factor 10 larger, and a plasma density of about 1.5 X 10^"^ 
cm~^. Since the plasma transits the scattering region in about 1 s, the scatterers must be 
replenished many times during the typical duration of an outburst of about a day. 

If, on the other hand, we assume that the outflow consists of a pair plasma then 
Lp lO"^" erg implies an outflow of 3 X 10"^^ pairs . If these pairs annihilate during 
the duration of the outburst, then the average annihilation flux at Earth is ~ 3 photons 

cm~^, which is much larger than the observed ~0.4 MeV line flux of (0.2-1) xlO~^ 
photons cm~^. Alternatively, if these pairs escape into the interstellar medium, the 
annihilation rate implied by the narrow 0.511 MeV line in the central Galactic bulge of 
~ 3 X 10"^^ (Ramaty, Skibo &: Lingenfelter 1994) limits the rate of occurance of day long 
outbursts to less then 1 in 1000 days which seems rather low. Finally, the minimum energy 
required to create these pairs is ~ 5 X lO"^" erg which exceeds the beamed gamma ray 
luminosity by a factor of 50. We thus believe that it is unlikely that the scattering medium 
is a pair plasma and consider only ordinary matter in the subsequent discussion. 

To produce line features as narrow as shown in Fig. 2, the temperature of the scat- 
tering medium must be less than about 5 keV. Since the line features are the result of 
electron recoil, the heating of the scattering medium is unavoidable. As a cooling mecha- 
nism we consider inverse Compton cooling by soft photons. We assume that the incident 
beam consists, in addition to gamma rays, also of UV photons with luminosity Luv The 
temperature is then kT ■qL-^mf.c^ / ALuv, where 77 is the fraction of L-^ deposited in the 
scattering medium. We find from our simulations that 77 0.2 for ry = 5 and so the 
temperature will be less than 5 keV if Luv ^5L-y. The line widths also depend on the 
coUimation of the incident beam. The results of Fig. 2 indicate that for spreads S/i in /i 
of 0.1 the lines have widths as seen in this figure. 
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We have estimated the 7-7 opacity due to the interaction of beamed photons with 
scattered photons. This opacity is smaller than that due to Compton scattering. We 
estimate that the pair production rate, dominated by 7-7 interactions, is about 10% of the 
angle integrated photon scattering rate, i.e. about lO"^"^ . Since the photon luminosity 
greatly exceeds the Eddington luminosity for pairs, we expect that most of these pairs will 
be blown out of the jets into the interstellar medium where they produce narrow 0.511 
MeV line emission. 

Mirabel &: Rodriguez (1994) have observed a brightening in the southern jet of 
1E1740. 7-2942 which they associated with bulk flow in the jet with apparent transverse 
velocity of 0.9 c. Our derived parameters, 7 1.2 and /x 0.45 yield an apparent velocity 
of only 0.65 c. Nonetheless, we believe that the radio observations do provide support for 
our model. The discrepancy in velocity could be due to several factors, including the time 
assumed by Mirabel &: Rodriguez (1994) at which the plasma blob was produced and the 
acceleration of the jet subsequent to the gamma ray line production. 

Although the geometry used here produces two lines, one from each jet, the model 
can accommodate single line spectra, for example that observed with SIGMA on October 
13-14, 1990 from 1E1740. 7-2942, if the beamed gamma ray emission or the jets themselves 
are asymmetric. We suggest that the spectrum observed with HEXAGONE (broad ~0.4 
and ^0.2 MeV features and a narrow 0.511 MeV line) represents the superposition of a 
double line spectrum from a compact object, possibly 1E1740. 7-2942, and narrow 0.511 
MeV line emission from the surrounding interstellar medium. 

We thank the referee, David Smith, for independently verifying our Monte-Carlo cal- 
culations and bringing to our attention an error in the original manuscript. One of us 
(CDD) acknowledges E. Churazov for pointing out that the Compton scattering of hard, 
beamed gamma rays leads to line features. 
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